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We present a new measurement of CP-violation parameters in B° —> Kg 7r°7 decay based on a 
sample of 275 x 10 6 BB pairs collected at the T(45 l ) resonance with the Belle detector at the KEKB 
energy-asymmetric e + e“ collider. One of the B mesons is fully reconstructed in the B° —> Kgn°-y 
decay. The flavor of the accompanying B meson is identified from its decay products. CP-violation 
parameters are obtained from the asymmetry in the distribution of the proper-time intervals between 
the two B decays. We obtain S K o w o 7 = —0.581^3® (stat) ± O.ll(syst) and v 4 fsr o„.o. r = +0.03 ± 
0.34(stat) ±0.11(syst), for the Kg-iv 0 invariant mass covering the full range up to 1.8GeV/c 2 . We 
also measure the CP-violation parameters for the case B° —> K *°(—> Kgn 0 )-y and obtain S K * o 7 = 
—0.791q fg(stat) ±0.10(syst) for A K *o^ fixed at 0. 

PACS numbers: 11.30.Er, 12.15.Hh, 13.25.Hw 


In the Standard Model (SM), CP violation arises 
from an irreducible phase, the Kobayashi-Maskawa (KM) 
phase Q, in the weak-interaction quark-mixing matrix. 
The phenomena of time-dependent CP violation in de¬ 
cays through radiative penguin processes such as b —> sy 
are sensitive to physics beyond the SM. Within the SM, 
the photon emitted from a B° (B°) meson is dominantly 
right-handed (left-handed). Therefore the polarization 
of the photon carries information on the original 6-flavor 
and the decay is, thus, almost flavor-specific. As a result, 
the SM predicts a small asymmetry .2, E] and any sig¬ 
nificant deviation from this expectation would be a man¬ 
ifestation of new physics. Recently it was pointed out 
that in decays of the type B° —► P°Q°7, where P° and 
Q° represent any CP eigenstate spin-0 neutral particles 
(e.g. P° = Kg and Q° = 7r°), many new physics effects 
on the mixing-induced CP violation do not depend on 
the resonant structure of the P°Q° system 0. In this 
Letter, we describe two measurements of CP asymme¬ 
tries: one with a Kgir 0 mass range restricted around the 
K*° [dlmass, and the other with an extended Kg n° mass 
range ja|. Whenever necessary, we distinguish these two 
by referring to them as K*° 7 and Kgir 0 ^, respectively; 
otherwise the analysis procedure is common to both mea¬ 
surements since it was first optimized for the K*° 7 and 
was extended for the Kgir 0 "/ later. 

At the KEKB energy-asymmetric e + e _ (3.5 on 
8.0GeV) collider 0, the T(45) is produced with a 


Lorentz boost of /3y = 0.425 along the z axis defined 
as antiparallel to the e + beam direction. In the decay 
chain T(45) —> B°B a —> / s i g /t ag , where one of the B 
mesons decays at time t s i g to a final state / s j g , which is 
our signal mode, and the other decays at time tt ag to a 
final state /tag that distinguishes between B° and B°, 
the decay rate has a time dependence given by 


V(A t) = 


g |At|/r s o 
4t b o 


1 + q 


S sin(AmdAf) 


+ Acos(AmdAt) 


(1) 


Here S and A are CP-violation parameters, tb 0 is the 
B° lifetime, A rrid is the mass difference between the two 
B° mass eigenstates, At is the time difference t s i g — f tag , 
and the 6-flavor charge q = +1 (—1) when the tagging 
B meson is a B° (B°). Since the B° and B° mesons are 
approximately at rest in the T (45 1 ) center-of-mass system 
(c.m.s.), At can be determined from the displacement in 
z between the / s i g and /t ag decay vertices: At ~ (z s i g — 
Zt ag )/(/37c) = A z/(/? 7 c). 

The Belle detector [g is a large-solid-angle magnetic 
spectrometer. A 2.0 cm radius beam pipe and a 3-layer 
silicon vertex detector (SVD1) were used for a 140/6 _1 
data sample containing 152 x 10 6 BB pairs, while a 1.5 cm 
radius beampipe, a 4-layer silicon detector (SVD2) 0 
and a small-cell inner drift chamber were used for an 
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additional 113/6 1 data sample that contains 123 x 10 6 
BB pairs for a total of 275 x 10 6 BB pairs. 

For high energy prompt photons, we select an isolated 
cluster in the electromagnetic calorimeter (ECL) that has 
no corresponding charged track, and has the largest en¬ 
ergy in the c.m.s. We require the shower shape to be 
consistent with that of a photon. In order to reduce the 
background from i r° or p, we exclude photons compati¬ 
ble with 7 r° —► 77 or 77 —■> 77 decays; we reject photon 
pairs that satisfy £ T o > 0.18 or £ v > 0.18, where /Vo 17-1 
is a 7T° (77) likelihood described in detail elsewhere [bQj. 
The polar angle of the photon direction in the labora¬ 
tory frame is restricted to the barrel region of the ECL 
(33° < 9 1 < 128°), but is extended to the end-cap re¬ 
gions (17° < 9 1 < 150°) for the second data sample due 
to the reduced material in front of the ECL. 

Neutral kaons (Kg) are reconstructed from two oppo¬ 
sitely charged pions that have an invariant mass within 
±6 MeV/c 2 (2cr) of the Kg nominal mass. The 7 r + 7 r _ ver¬ 
tex is required to be displaced from the interaction point 
(IP) in the direction of the pion pair momentum El- 
Neutral pions (77°) are formed from two photons with 
the invariant mass within ±16 MeV/c 2 (3cr) of the 7 r° 
mass. The photon momenta are then recalculated with 
a 7T° mass constraint and we require the momentum of 
7T° candidates in the laboratory frame to be greater than 
0.3 GeV/c. The Kgir 0 invariant mass, M K o^o, is required 
to be less than 1.8 GeV/c 2 . 

B° mesons are reconstructed by combining Kg, 7 r° and 
7 candidates. We form two kinematic variables: the en¬ 
ergy difference A E = Eg ms - — and the beam- 

energy constrained mass M bc = y/(Ebeam') 2 — (Ps m ' s ') 2 ) 
where -E^eam’ is the beam energy, Eg ms - and p c g ms ' are 
the energy and the momentum of the candidate in the 
c.m.s. Candidates are accepted if they have M bc > 
5.2 GeV/c 2 and \ AE\ < 0.5 GeV. 

We reconstruct B + —■» Kgir + "/ candidates in a sim¬ 
ilar way as the B° —> Kgir 0 "/ decay in order to re¬ 
duce the cross-feed background from B + —> Kgir + "/ in 
B° —> Kgir 0 "/. The B + —> K°ir + "/ events are also used 
for various cross-checks. For a 7 r + candidate, we require 
that the track originates from the IP and that the trans¬ 
verse momentum is greater than 0.1 GeV/c. We also re¬ 
quire that the 7 r + candidate cannot be identified as any 
other particle species ( K + ,p + ,e + , and p + ). 

Candidate B + —> Kgir + "/ and B° —► Kgir 0 "/ decays are 
selected simultaneously; we allow only one candidate for 
each event. The best candidate selection is based on the 
event likelihood ratio 77. s / b that is obtained by combining 
a Fisher discriminant T , which uses the extended 
modified Fox-Wolfram moments El as discriminating 
variables, and cos9h defined as the angle between the 
B meson momentum and the daughter Kg momentum 
in the rest frame of the Kgir system. We select the can¬ 
didate that has the largest 77. s / b . For the exclusive B° —> 


K *°7 study, we further require 0.8 GeV/c 2 < M K o n o < 
1.0 GeV/c 2 after the best candidate selection. The sig¬ 
nal region is defined as —0.2 GeV < A E < 0.1 GeV and 
5.27 GeV/c 2 < M bc < 5.29 GeV/c 2 . 

We use events outside the signal region as well as large 
Monte Carlo (MC) samples to study the background 
components. The dominant background is from con¬ 
tinuum light quark pair production (e + e~ —*► qq with 
q — u,d, s, c), which we refer to as qq hereafter. We re¬ 
quire 77. s /b > 0.5 to reduce the qq background; after ap¬ 
plying all other selection criteria, this rejects 79% of the 
qq background while retaining 88 % of the signal. Back¬ 
ground contributions from B decays, which are consider¬ 
ably smaller than qq, are dominated by cross-feed from 
other radiative B decays including B + —► K°ir + "/. 

The 6 -flavor of the accompanying B meson is identified 
from inclusive properties of particles that are not associ¬ 
ated with the reconstructed signal decay. The a lgor ithm 
for flavor tagging is described in detail elsewhere [14J. We 
use two parameters, q defined in Eq. CD and r, to rep¬ 
resent the tagging information. The parameter r is an 
event-by-event flavor-tagging dilution factor that ranges 
from 0 to 1 ; r = 0 when there is no flavor discrimination 
and r = 1 implies unambiguous flavor assignment. It is 
determined by using MC data and is only used to sort 
data into six r intervals. The wrong tag fraction w and 
the difference Aw in w for the B° and B° decays are 
determined for each of the six r intervals from data El 

The vertex position of the signal-side decay is recon¬ 
structed from the Kg trajectory with a constraint on the 
IP; the IP profile (a x — 100 pm, a y ~ 5 pm, 07 ~ 3 mm) 
is convolved with the finite B flight length in the plane 
perpendicular to the z axis. Both pions from the K° 
decay are required to have enough hits in the SVD in 
order to reconstruct the Kg trajectory with high resolu¬ 
tion. The reconstruction efficiency depends not only on 
the Kg momentum but also on the SVD geometry. The 
efficiency with SVD2 (55%) is higher than with SVD1 
(41%) because of the larger detection volume. The other 
(tag-side) B vertex determination is the same as that for 
the B° —> (j)K° s analysis El 

Figure Q] shows the M bc (A E) distribution for the re¬ 
constructed Kgir 0 "/ candidates within the A E (M bc ) sig¬ 
nal region after flavor tagging and vertex reconstruction. 
The signal yield is determined from an unbinned two- 
dimensional maximum-likelihood fit to the AE-M\, C dis¬ 
tribution. The fit region for the B° —> K *°7 analysis 
is 5.20 GeV/c 2 < M bc < 5.29 GeV/c 2 and -0.5 GeV < 
A E < 0.5 GeV. For the B° -> K° s ir 0 "/ analysis, the A E 
fit region is narrowed to —0.4 GeV < A E < 0.5 GeV in 
order to avoid other B background events that populate 
the low A E region. The signal distribution is represented 
by a smoothed histogram obtained from MC simulation 
that accounts for a small correlation between M bc and 
A E. The background from B decays is also modeled with 
a smoothed histogram obtained from MC events. For the 
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FIG. 1: (a) Mbc distributions within the A E signal region 
and (b) A E distributions within the Mbc signal region for 
B° ->■ 7r°7. Solid curves show the fit to signal plus back¬ 

ground distributions. Lower (upper) dashed curves show the 
background contributions from qq (qq and B decays). 
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FIG. 2: M K o n o distribution of the B° —> K%-k 0 ^ signal yield 
obtained from the fit to the AE-Mhc distribution. Back¬ 
ground shape parameters are fixed at the values obtained with 
the whole M K o n o region. 


qq background, we use the ARGUS parameterization m 
for Mbc and a second-order polynomial for A E. Normal¬ 
izations of the signal and background distributions and 
the qq background shape are the five free parameters in 
the fit. We observe 501 (215) Kgir 0 ^ (K* 0r y) candidates 
in the signal box, which decrease to 227 (92) after the 
flavor tagging and the B vertex reconstruction, and ob¬ 
tain 105 ± 14 (57 ± 9) signal events by the fit. Figure 0 
shows the signal yields in different M K o n o regions. For 
each bin, a fit to the AE-M^ C distribution is performed, 
and the obtained signal yield is plotted. 

We determine S and A from an unbinned maximum- 
likelihood fit to the observed At distribution. The proba¬ 
bility density function (PDF) expected for the signal dis¬ 
tribution, V s i g (At',S,A,q,w,Aw), is given by the time 
dependent decay rate [Eq. (0)] modified to incorporate 
the effect of incorrect flavor assignment. The distribu¬ 
tion is convolved with the proper-time interval resolution 
function R s i g , which takes into account the finite vertex 
resolution. The parametrization of R s is the same as 
that used for the B° —> Kgn 0 decay ll|. R s i g is first 
derived from flavor-specific B decays Il6l and modified 
by additional parameters which rescale vertex errors to 
account for the fact that there is no track directly origi¬ 
nating from the IP. 

For each event, the following likelihood function is eval¬ 


uated: 
Pi = 


/ +oo 

-oo 


/sig'Psig ( At') f?sig ( At, - At') 


± (1 - /sig)'Pbkg(Af , )i?bkg(Ati - At') 
± /ol-Pol(Afi), 


d(At') 


( 2 ) 


where P 0 i is a Gaussian function that represents a small 
outlier component with fraction f Q i E2 The signal prob¬ 
ability / s ig is calculated on an event-by-event basis from 
the function which we obtained as the result of the two- 
dimensional AU-Mbc fit for the signal yield extraction. 
A PDF for background events, Pbkg, is modeled as a sum 
of exponential and prompt components, and is convolved 
with a sum of two Gaussians which represent the resolu¬ 
tion function i?bkg for the background. All parameters 
in Pbkg and I?bkg are determined by a fit to the At dis¬ 
tribution of a background-enhanced control sample, i.e. 
events outside of the AE-M\> C signal region. We fix t b o 
and Arrid at their world-average values [13 ■ 

The only free parameters in the final fit for the B° —* 
Kgir°'y decay are S K o 7 r o 7 and 4 ^- 0 ^° 7 , which are de¬ 
termined by maximizing the likelihood function L = 
I], Pi{Ati\ S, A) where the product is over all events. For 
the B° —► K *°7 decay, we fix Ak* o 7 at zero and perform 
a single parameter fit for Sk*°^- This is based on the pre¬ 
vious measurement of the partial rate asymmetry in 
B —> K *7 decays, which yielded Ak* o 7 consistent with 
zero. The error of the previous measurement of Ak *° 7 is 
included in the systematic error of Sk* o 7 . We obtain 

Sk*o-j = —0.79(toJo(stat) ± O.lO(syst), 
o 7 = -0.58t^|(stat) ± O.ll(syst), 

7 = +0.03 ± 0.34(stat) ± O.ll(syst). 

We define the raw asymmetry in each At bin by 
(A(j = + 1 - fVq=-i)/(Af g =+i + N q= ^ 1 ), where AT g=+1( _ 1 ) 
is the number of observed candidates with q = + 1 (— 1 ). 
Figure 0 shows the raw asymmetries for the KgTT 0r y and 
K* 0r y decays in two regions of the flavor-tagging param¬ 
eter r. Note that these are simple projections onto the 
At axis, and do not reflect other event-by-event informa¬ 
tion (such as the signal fraction, the wrong tag fraction 
and the vertex resolution), which is in fact used in the 
unbinned maximum-likelihood fit for S and A. 

Primary sources of the systematic error are (1) un¬ 
certainties in the resolution function (±0.05 for S K g w o 7 , 
±0.04 for A k o t o 7 and ±0.05 for Sk* o 7 ), ( 2 ) uncertainties 
in the vertex reconstruction (±0.05 for S^o^o^, ±0.06 
for Ak°-k o 7 and ±0.06 for <S/f*o 7 ), (3) uncertainties in 
the background fraction (±0.06 for 5^o w o 7 , ±0.04 for 
A^-o^o 7 and ±0.02 for Sk* o 7 ), and (4) uncertainties in 
the background At distribution (±0.05 for l S^-o T o 7 , ±0.04 
for 4 if o I o 7 and ±0.02 for Sk *° 7 )- Effects of tag-side in¬ 
terference JliJ contribute ± 0.01 for 5 A -o^o 7 , ±0.06 for 
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FIG. 3: Asymmetry in each At bin with 0 < r < 0.5 (top) 
and with 0.5 < r < 1.0 (bottom) for (a) B° — > K*° 7 and for 
(b) B° — > Kgn 0 ^. Solid curves show the results of the un¬ 
binned maximum-likelihood fits. The dashed curves show the 
null asymmetry case (S = A = 0). The dilution of a possible 
asymmetry is smaller for the 0.5 < r < 1.0 subsample, be¬ 
cause of the lower probability of incorrect flavor assignment. 

A k o 7r o 7 and ±0.01 for Sk *° t - Also included are effects 
from uncertainties in the wrong tag fraction and physics 
parameters (A rrid, t b 0 and Ak* t> 7 ). A possible bias in 
the fit is estimated by using MC data. The total sys¬ 
tematic error is obtained by adding these contributions 
in quadrature. 

Various cross-checks of the measurement are per¬ 
formed. We apply the same fit procedure to the B° —> 
J/ipKg sample without using J/ip daughter tracks for 
the vertex reconstruction. We obtain Sj/^ K ° s = ±0.68 ± 
O.lO(stat) and Aj/^k° = ±0.02 ± 0.04(stat), which are 
in good agreement with the world average values 0 - 
We perform a fit to B + —> Kgir+'y (B + —> K * + 7 —> 
A'Ott+ 7 ), which is a good control sample of the B° —> 
A'^ 7 t °7 (B° —*■ K *°7 —► decay, without using 

the primary 7 r + for the vertex reconstruction. The result 
is consistent with no CP asymmetry, as expected. Life¬ 
time measurements are also performed for these modes, 
and values consistent with the world-average are ob¬ 
tained. Ensemble tests are carried out with MC pseudo¬ 
experiments; we find that the statistical errors obtained 
in our measurements are all within the expectations from 
the ensemble tests. A fit to the B° —► KgTT 0 j sub¬ 
sample that excludes the K*° mass region yields S = 
—0.39lo 52 ( s tat) and A = ±0.10 ± 0.51(stat). The result 
is consistent with the results for B° —> AT *°7 and for the 
full B° —> ATg 7 r °7 sample. 

In summary, we have performed a new measurement of 
the time-dependent CP asymmetry in the decay B° —> 
ATg7r°7 based on a sample of 275 x 10 6 BB pairs. Two 
regions of the K B n° invariant mass are used: between 
0.8GeV/c 2 and 1.0GeV/c 2 for the K*° resonance, and 
the full region up to 1.8GeV/c 2 including other res¬ 
onances and non-resonant contributions. We obtain 
CP-violation parameters in the B° —> Kg 7r°7 decay 
S K o■ n ° 1 = -0-58to:3s( stat ) ± O.ll(syst) and A K o 7r o 7 = 
±0.03 ± 0.34(stat) ± O.ll(syst), and in B° —> K*°(—> 


A'g7 r °)7 decay 5/y.o 7 = —0.79lo]5o( s tat) = tO-10(syst) with 
Ak *° 7 fixed at zero. The value of iSr-»o 7 is consistent with 
the measurement by BABAR [ill ] . CP-violation parame¬ 
ters in the three-body B° —> Kg7r 0 "f decay are measured 
for the first time. The two results are consistent with 
each other and with no CP asymmetry, while the sta¬ 
tistical error is considerably reduced by the inclusion of 
these additional events in the B° — > Kgir 0 ^ analysis. 

We thank the KEKB group for the excellent operation 
of the accelerator, the KEK cryogenics group for the ef¬ 
ficient operation of the solenoid, and the KEK computer 
group and the Nil for valuable computing and Super- 
SINET network support. We acknowledge support from 
MEXT and JSPS (Japan); ARC and DEST (Australia); 
NSFC (contract No. 10175071, China); DST (India); the 
BK21 program of MOEHRD and the CHEP SRC pro¬ 
gram of KOSEF (Korea); KBN (contract No. 2P03B 
01324, Poland); MIST (Russia); MHEST (Slovenia); 
SNSF (Switzerland); NSC and MOE (Taiwan); and DOE 
(USA). 


* on leave from Nova Gorica Polytechnic, Nova Gorica 

[1] M. Kobayashi and T. Maskawa, Prog. Theor. Phys. 49, 
652 (1973). 

[2] D. Atwood, M. Gronau and A. Soni, Phys. Rev. Lett. 79, 
185 (1997). 

[3] B. Grinstein, Y. Grossman, Z. Ligeti and D. Pirjol, Phys. 
Rev. D 71, 011504 (2005). 

[4] D. Atwood, T. Gershon, M. Hazumi and A. Soni, Phys. 
Rev. D 71, 076003 (2005). 

[5] Throughout this paper, K* denotes K*(892). The inclu¬ 
sion of the charge conjugate decay mode is implied. 

[ 6 ] Possible strong phase dependence within the SM was 
pointed out by Grinstein et al. .jj|. In the case CP vi¬ 
olation depends on the K^tv 0 mass, our measurement for 
Kg 7 r °7 gives an efficiency-weighted average of CP asym¬ 
metries over the extended mass range. 

[7] S. Kurokawa and E. Kikutani, Nucl. Instr. and Meth. A 
499, 1 (2003). 

[ 8 ] Belle Collaboration, A. Abashian et al., Nucl. Instr. and 
Meth. A 479, 117 (2002). 

[9] Y. Ushiroda (Belle SVD2 Group), Nucl. Instr. and Meth. 
A 511, 6 (2003). 

[10] Belle Collaboration, P. Koppenburg et al., Phys. Rev. 
Lett. 93, 061803 (2004). 

[11] Belle Collaboration, K. Abe et al, hep-ex/0409049 

[12] R. A. Fisher, Annals Eugen. 7, 179 (1936). 

[13] Belle Collaboration, K. Abe et al, Phys. Rev. Lett. 91, 
261801 (2003). 

[14] H. Kakuno, K. Hara et al, Nucl. Instr. and Meth. A 533, 
516 (2004). 

[15] ARGUS Collaboration, H. Albrecht et al., Phys. Lett. B 
241, 278 (1990). 

[16] Belle Collaboration, K. Abe et al., hep-ex/0408111 

[17] Belle Collaboration, K. Abe et al., Phys. Rev. D 71, 
072003 (2005) ; H. Tajima et al, Nucl. Instrum. Methods 
Phys. Res., Sect. A 533, 370 (2004). 

[18] S. Eidelman et al., Phys. Lett. B 592, 1 (2004). 


















6 


[19] Belle Collaboration, M. Nakao et al, Pliys. Rev. D 69, [21] BABAR Collaboration, B. Aubert et al, Phys. Rev. Lett. 

112001 (2004). 93, 201801 (2004). 

[20] O. Long, M. Baak, R. N. Cahn and D. Kirkby, Phys. 

Rev. D 68, 034010 (2003). 


